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The  Speckle  Holography  technique  uses  a  series  of  short  exposure  images  and  associated  wave  front  sensor 
measurements  to  reconstruct  astronomical  images.  The  reconstruction  technique  effectively  boosts  the  amplitude 
of  high  spatial  frequencies  to  allow  for  finer  resolution  of  astronomical  objects.  A  simulation  package  that 
calculates  the  optical  transfer  function  (OTF)  of  a  telescopic  system  was  previously  developed  for  the  Phillips 
Laboratory.  This  software  was  modified  slightly  to  provide  the  Speckle  Holography  system  transfer  function 
(STF)  for  a  telescopic  system  under  varying  seeing  conditions  as  specified  by  the  ratio  of  the  subaperture 
separation  {La,)  to  the  Fried  parameter  (r,),  and  by  the  photon  count  per  subaperture  per  short  exposure 
image.  The  thesis  analyzes  two  different  approaches  to  the  STF.  In  addition  to  generating  STF  curves  for 
varying  seeing  conditions,  the  signal-to-noise  ratio  (SNR)  for  a  single  sampling  of  short  exposure  imap  and  wave 
front  measurements  is  calculated  according  to  the  formula  SNR  =  STF/-^var,  where  var  is  the  variance  of  the 
STF.  The  thesis  results  indicate  the  STF  ic  severely  degraded  by  a  bias  that  originates  from  the  estimate  of 
the  OTF,  which  is  used  in  the  STF  calculation.  The  problem  is  found  to  be  increasingly  severe  for  improved 
wave  front  sensor  sampling  (denser  wave  front  measurements)  and  low  photon  counts.  An  attempt  is  made  to 
characterize  non-biased  regions  of  the  STF. 
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Tlu'  Spccklo  llolograpliy  ttThiii(|ii(‘  usi's  a  writ's  of  short  ('xposiiri'  iinaps  ainl  associali'd  wavi- 
front  st’nsor  nirasuri’iin-nts  to  rt'constriict  astronomiral  itnagi-s.  I'lii’  roronstriirtion  ti'chniipii'  I'ffcc- 
tivoly  boosts  tho  aiiiplitiidt' of  high  spatial  frt'<iii<'iiri<’s  to  allow  for  (im-r  ri’sohitioii  of  aslronoiniral 
ohjorts.  A  siinnlation  parkago  that  ralrulatfs  tho  optical  transfor  fiiiirtioii  (O  TF)  of  a  ti-lcsct'pic 
systrni  was  previously  dowlopod  for  the  Phillips  Laboratory.  This  software  was  nioclilied  slightly 
to  provide  the  Speckle  Holography  system  transfer  function  (STF)  for  a  telescopic  system  under 
varying  sc'eing  conditions  as  specified  by  the  ratio  of  the  subapi'rtiire  separation  (/-,i. )  to  the  Fried 
parameti-r  (r,,).  and  by  the  photon  count  per  siibapertiire  per  short  exposun*  image.  The  thesis 
aualyz('s  two  different  approaches  to  the  STF.  In  addition  to  generating  STF  ciirvt's  for  varying 
sei'ing  conditions,  the  signal-to-noi.se  ratio  (SNR)  for  a  single  sampling  of  short  exposure  image  and 
wave  front  measurements  is  calculated  according  to  the  formida: 


SNRr: 


STF 

\/var’ 


whore  var  is  the  variance  of  the  STF.  The  the.sis  results  indicate  the  STF  is  .wverely  degraded 
by  a  bias  that  originatt's  from  the  estimate  of  the  OTF,  which  is  uwd  in  the  S'FF  calculation. 
The  problem  is  found  to  be  increasingly  severe  for  improved  wave  front  wn.sor  sampling  (denwr 
wave  front  measurements)  and  low  photon  counts.  An  attempt  is  made  to  characterize  non-biased 
regions  of  the  STF. 


IX 


SELF-REFEHKN’rKI)  SPECKLE  HOLOGRAPHY  (SRSH) 
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/.  Aiviosphcrir  Turbulence  and  Ihe  SliSH  Technique, 
t.i  Introduelton. 

Prior  to  the  space-hased  Huhhle  telescope,  grou;  J  based  telescopes  were  the  only  means  for 
ohsi'rving  distant  stars.  Although  a  space  bas<'d  bdescopi-  has  the  advantage  of  not  having  to  “look” 
through  the  atmosphere,  maintenance  costs  and  si/<>  limitations  for  spaced  based  systems  ensure 
that  groimd-ba.sed  telescopes  will  continue  to  remain  a  vital  tool  for  exploring  the  universe. 

From  a  diffraction  theory  perspective,  a  ground  based  telescope  should  be  built  as  large  as 
possible  since  the  larger  the  leh'scope  diameter,  the  finer  the  angular  resolution  of  the  telescope  and 
the  smaller  the  object  the  telescope  can  image.  It  turns  out  that  the  resolution  of  large  diameter 
telescopes  is  limited  by  the  effects  of  atmospheric  turbulence.  The  turbulence  induces  random 
spatial  and  temporal  fluctuations  on  wave  fronts  prior  to  imaging  by  the  telescope,  blurring  the 
image  and  decreasing  the  angular  re.solution  from  the  theoretical  limit. 

Wave  front  measurements  and  signal  processing  can  be  used  to  partially  overcome  the  effects 
of  turbulence  on  images.  In  this  thesis,  a  FORTRAN  computer  simulation  is  used  to  examine 
the  imaging  performance  of  the  self- referenced  speckle  holography  (SRSH)  image  reconstruction 
technique.  This  technique  shows  great  promise,  with  imaging  results  approaching  the  theoretical 
angular  resolution  limit  of  large  diameter  tele.scopes  to  which  it  is  applied  [1:1].  In  this  chapter, 
the  reader  is  introduced  to  the  problem  of  imaging  through  atmq.spheric  turbulence,  followed  by  an 
introduction  to  the  SRSH  technique,  and  general  thesis  objectivt's.  In  Chapter  II,  the  general  theory 
of  the  SRSH  technique  is  p  esented  along  with  assumptions  and  defining  equations  for  the  thesis. 
The  FORTRAN  computer  simulation  is  described  in  some  detail  in  Chapter  III.  The  results  for 
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this  thesis  effort  are  presented  in  Chapter  IV  and  conclusions  and  reccmniendations  are  presented 
in  Chapter  V. 

1.2  Almospheric  Turbulence. 

1.2.1  Turbulent  eddies.  Atmospheric  turbulence  is  generated  by  uneven  heating  and 
cooling  of  the  earth’s  atmosphere.  This  process  results  in  randomly  generated  pockets  of  turbulence 
or  eddies,  each  of  wiiich  has  an  index  of  refraction  defined  by  [2:56}: 


where  n  is  the  index  of  refraction,  c  (m/s)  is  the  speed  of  light  in  a  vacuum,  and  v  (m/s)  is  the 
velocity  of  propagation  of  light  in  the  medium.  Each  pocket  of  turbulence  can  have  a  different 
index  of  refraction  from  other  pockets  of  turbulence.  These  pockets  occur  randomly  in  a  temporal 
and  spatial  sense,  and  so  induce  random  temporal  and  spatial  fluctuations  on  wave  fronts  passing 
through  the  turbulent  region,  causing  image  distortions  [3:4527]. 

Figure  1.1  shows  a  wave  front  prior  to  and  after  encountering  turbulence  while  enroute  to  a 
ground  based  telescope.  Light  Lorn  a  star  is  es,sentially  planar  by  the  time  it  reaches  the  earth’s 
atmosphere  due  to  the  distance  of  propagation.  In  traveling  through  the  atmosphere  to  a  ground 
based  telescope,  the  planar  wave  front  encounters  the  randomly  distributed  turbulent  eddies.  The 
overall  delay  experienced  by  a  specific  point  on  the  wave  front  will  depend  on  which  turbulent 
eddies  that  point  encounters  as  it  passes  through  the  turbulent  atmosphere. 

1.2.2  Minimum  angular  resolution  in  turbulent  conditions.  The  net  effect  of  this  tur¬ 
bulence  is  to  reduce  the  angular  resolution  of  a  large  telescope.  In  theory,  the  diffraction  limited 
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Figure  1.1.  Wave  front  degradation  due  to  atmospheric  turbulence, 
angular  resolution  of  a  telescope  is  given  by  [2:422]; 

^<»min=122^.  (1.2) 

where  's  Ih®  angular  resolution  of  the  telescope,  A  (m)  is  the  wavelength  of  the  incident 

light,  and  D  (m)  is  the  telescope  diameter.  Equation  (1.2)  implies  the  angular  resolution  imp'oves 
with  increasing  telescope  diameter. 

Fried  derived  a  parameter  in  1966  to  describe  the  point  at  which  increasing  atmospheric 
distortion  limits  the  resolution  performance  of  a  telescope  [4:33].  The  so-called  Fried  coherence 
diameter  or  Fried  parameter,  Vg,  is  a  measure  of  the  lateral  spatial  coherence  properties  of  incident 
wave  fronts  [5:10-13].  The  Fried  parameter  is  roughly  on  the  order  of  the  size  of  the  sloped  regions 
induced  on  a  wave  front.  The  Fried  parameter  is  a  function  of  the  wavelength  of  the  light,  the 
propagation  path  length  through  the  turbulence,  the  zenith  angle  (angle  with  the  vertical),  and  the 
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strength  of  the  turbulence.  If  the  telescope  diameter  is  less  than  the  Fried  parameter,  the  telescope 
operates  with  essentially  diffraction  limited  performance.  When  the  telescope  diameter  exceeds  the 
Fried  parameter  (as  is  generally  the  case  for  large  diameter  telescopes  at  visible  wavelengths),  the 
angular  resolution  achievable  is  aiiproximated  by  replacing  D  in  equation  (1.2)  by  Va  to  give  [6:14]: 

=  I.22A.  (1.3) 

The  minimum  angular  resolution  of  large  diameter  ground  based  telescopes  is  generally  limited  to 
around  1  arcsecond  at  visible  wavelengths  [7:1-2].  For  example,  the  2.4  meter  diameter  Hubble 
telescope  has  a  diffraction  limit  of  0.03  arcsecond  in  space.  On  the  ground,  this  same  telescope 
would  be  limited  to  angular  resolutions  on  the  order  of  0.5  to  1.0  arcseconds  [6:1,14],  which  is  16 
to  30  times  worse  than  the  theoretical  limit. 

The  Fried  parameter,  To,  typically  ranges  from  3  centimeters  for  poor  seeing  conditions  to  20 
or  30  centimeters  for  good  seeing  conditions.  Its  value  is  site  dependent,  since  the  altitude  of  an 
observatory  site  determines  how  much  of  the  atmosphere  the  light  propagates  through.  Its  value  is 
also  weather  dependent  [6:14].  Typically,  atmospheric  turbulence  limits  large  diameter  telescopes  to 
the  same  angular  resolution  as  smaller  telescopes  on  the  order  of  10  to  20  centimeters  in  diameter. 

1.3  Historical  Development  of  the  SRSH  Technique. 

Several  image  recovery  techniques  have  been  developed  that  involve  computer  post  processing 
of  a  series  of  short  exposure  images  to  recover  an  improved  image  of  the  original  object.  The 
reason  for  using  short  exposure  images  is,  according  to  Goodman,  to  “‘freeze’  the  atmospheric 
degradations,  thus  eliminating  any  time  averaging  effects  [8:402].”  A  long  exposure  image  is  blurred 
by  the  changing  atmospheric  conditions  causing  the  image  to  move  around  during  the  exposure. 
These  image  changes  occur  on  the  order  of  milliseconds  [6].  An  example  of  a  long  exposure  image 
is  shown  in  Figure  1.2  (a).  A  short  exposure  image,  see  Figure  1.2  (b),  contains  diffraction  limited 
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spots  or  speckles.  These  speckles  result  “from  interference  of  a  signal  with  itself  or  with  part  of 
the  actual  imaging  or  optical  system  [1:1].”  A  series  of  short  exposure  images  can  be  processed 
to  reconstruct  a  near  diffraction  limited  image  as  evidenced  by  Labeyrie’s  Speckle  Interferometry 
technique,  the  Knox-Thompson  algorithm  technique,  and  the  Triple  Correlation  technique.  All  of 
these  image  post  processing  techniques  can  effectively  reconstruct  images  out  to  near  the  diffraction 
limit  of  the  telescope.  The  problem  with  these  techniques  is  they  can  require  thousands  of  short 
exposure  images  [9:64]. 

An  alternative  to  post  processing  is  what  Hardy  calls  active  optics  or  adaptive  optics  com¬ 
pensation  [10:654].  For  this  technique,  distorted  wave  fronts  are  compensated  real  time  by  use  of 
a  controllable  corrective  surface  whose  refractive  or  reflective  properties  are  changed  to  cancel  the 
incident  wave  front  distortion.  The  technique  consists  of  measuring  the  phase  distortion  along  the 
wave  front  and  then  applying  the  appropriate  controls  to  the  refractive  or  reflective  surface.  The 
three  parts  to  such  a  system  are  the  controllable  corrective  surface,  wave  front  sensor  device,  and 
information  processing  device  required  to  generate  control  signals  for  the  corrective  surface. 

One  type  of  commonly  used  wave  front  .sensor  device  is  the  Shack-Hartmann  wave  front 
sensor  (WFS),  which  h2is  the  ideal  capability  of  detecting  local  wave  front  slopes  at  low  light  levels 
[1:11].  Figure  1.3  shows  an  array  of  subapertures  in  a  Shack-Hartman  WFS.  The  telescope  aperture 
of  length  D  is  subdivided  into  smaller  apertures  of  dimension  L.  Each  of  these  small  apertures, 
or  subapertures,  has  a  lens  that  focuses  that  portion  of  the  wave  incident  on  the  subaperture 
onto  a  detector.  The  local  slope  of  the  incident  wave  front  is  sensed  separately  in  the  x  and  y 
directions  by  measuring  the  distance  the  light  is  focused  off  axis  as  shown  in  Figure  1.4.  The  Shack- 
Hartman  WFS  provides  the  key  capability  of  local  tilt  measurements  for  the  SRSH  technique,  and 
its  development  resulted  in  practical  application  of  the  SRSH  technique  in  image  reconstruction. 

Fontanella first  proposed  enhancing  post  processing  image  reconstruction  performance  through 
a  new  technique  that  used  simultaneous  monochromatic  short  exposure  speckled  images  and  wave 
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Figure  1.2.  Examples  of  (a)  long  and  (b)  short  exposure  images.  [8:403] 
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Figure  1.3.  Shack-Hartman  wave  front  sensor  (VVFS)  subaperture  configuration. 

front  measurements  using  the  Shack-Hartman  sensor  [9:64].  This  new  hybrid  approach  combined 
short  exposure  image  post  processing  with  the  adaptive  optics  WFS  measurements  to  remove  the 
effect  of  the  turbulent  atmosphere.  This  technique  is  called  Speckle  Holography  due  to  the  use  of 
speckled  short  exposure  images  such  as  that  shown  in  Figure  1.2  (b),  and  the  use  of  a  reference  wave 
front  (similar  to  holography  recordings)  from  a  point  source  for  phase  measurements  [11:4527].  The 
SRSH  technique  differs  only  from  Speckle  Holoff'aphy  in  that  for  the  SRSH  technique,  the  object 
itself  acts  m  its  own  reference  for  phase  mecisurements. 

1.4  General  Thesis  Objectives. 

This  thesis  will  address  the  performance  of  the  SRSH  technique  by  first  investigating  the 
validity  of  results  from  earlier  work  by  VonNiederhausern  that  characterized  the  ensemble  average 
Fourier  domain  system  transfer  function  (STF)  for  the  SRSH  technique  [16].  The  STF  is  the  Fourier 
domain  trans^'er  function  for  the  inversion  from  the  degraded  image  spectrum  to  the  estimated 
object  spectrum.  The  SRSH  technique  removes  turbulence  degradations  in  the  image  by  filtering  in 


Figure  1.4.  Side  view  of  Shack-Hartnian  wave  front  sensor  (WFS)  lens  and  detector  array. 


the  Fourier  domain.  The  STF  results  from  this  earlier  work  were  characterized  for  variations  in  the 
Fried  parameter  (r#)  and  the  strength  of  the  incident  light  (average  photon  count  per  subaperture 
per  short  exposure  image).  This  thesis  uses  a  FORTRAN  simulation  to  generate  STF  results  by 
performing  calculations  on  a  limited  numbers  of  short  exposure  images  and  WFS  measurements. 
This  thesis  will  consider  two  different  approaches  to  calculating  the  STF  as  defined  in  Chapter  II. 
Additionally,  this  thesis  presents  results  for  a  signal  to  noise  ratio  (SNR)  analysis  that  characterizes 
the  single  frame  STF  imaging  performance  as  a  function  of  spatial  frequency.  An  unexpected  STF 
bias  problem  is  also  investigated  and  characterized  for  this  thesis. 


A  discussion  of  the  hardware  required,  the  basic  theory  of  the  SRSH  technique,  and  assump¬ 


tions  and  defining  equations  for  the  thesis  appears  in  the  next  chapter. 


1-8 


II.  General  Theory  of  the  SRSII  Technique. 


2.1  Introduction. 

This  thesis  presents  results  from  a  FORTRAN  computer  simulation  of  a  telescope  configured 
with  the  hardware  shown  in  Figure  2.1.  As  background  for  the  SRSH  simulation,  the  basic  Fourier 
theory  for  a  simple  lens  equivalent  model  is  presented.  This  simple  system  configuration  is  basicly 
the  optical  system  used  by  the  simulation.  Next,  the  defining  equation  of  the  SRSH  technique,  the 
sjstem  transfer  function  (STF),  is  derived  using  the  simple  lens  model.  Finally,  the  a-ssumptions 
and  defining  equations  for  the  thesis  are  presented. 

2.2  Hardware  Configuration  for  the  SRSH  Technique. 

As  shown  in  Figure  2.1,  the  SRSH  technique  requires  a  fast  shutter  to  control  exposure  times, 
a  beam  splitter  to  split  an  incident  wave  front  for  simultaneous  short  e.xposure  image  and  wave 
front  measurements,  a  wave  front  sensor  for  wave  front  measurements  (Shack-Hartman  sensor),  a 
filter  to  pass  quasimonochromatic  light,  a  camera  for  recording  the  short  exposure  images,  and 
a  computer  to  construct  an  estimate  of  the  original  object  using  wave  front  measurement  and 
short  exposure  image  data  [11:4528].  The  image  is  a  two  dimensional  intensity  distribution,  i(x,  y), 
while  the  wave  front  measurements  are  actually  slope  measurements  that  are  used  to  construct  an 
estimate,  0(x,j/),  of  the  phase  of  the  incident  wave  front  phase,  <j>{i:,y).  The  simulation  uses  an 
adaptive  optics  control  law  to  construct  the  phase  estimate,  y).  An  estimate  of  the  object  being 
imaged,  d{x,y),  is  obtained  by  applying  the  SRSH  algorithm  over  many  frames  of  simultaneous 
phase  estimates  and  short  exposure  images. 

2..S  Lens  Equivalent  Model. 

For  most  astronomical  scenarios  involving  large  diameter  telescopes,  and  good  "seeing”  con¬ 
ditions  where  clouds  and  other  aerosols  are  not  present,  the  atmospheric  distortion  induces  mainly 
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Figure  2.1.  Typical  SRSH  hardware  configuration. 


phase  distortions  [8:362].  This  result  occurs  when  the  turbulence  through  which  the  telescope  is 
viewing  the  object  is  constrained  to  what  is  called  the  near  field.  For  zenith  angles  close  to  zero  and 
a  large  diameter  telescope,  one  can  expect  less  than  a  4  percent  error  in  calculations  due  to  neglect¬ 
ing  the  amplitude  distortions  [12:588-589].  The  near  field  assumption  holds  for  the  FORTRAN 
simulation  since  the  telescope  diameter  is  sufficiently  large  and  the  zenith  angle  is  zero.  So,  the 
turbulent  atmosphere  for  short  exposure  image  formation  can  be  modeled  by  a  single  phase  screen 
that  introduces  delays  across  the  incident  wave  front  [11:4527].  A  sample  wave  front  coming  from 
such  a  phase  screen  i])  shown  in  the  greyscale  plot  (81  by  81  pixels)  in  Figure  2.2  with  bright  pi.tels 
corresponding  to  small  delays  and  dark  pixels  to  large  delays.  Figure  2.3  shows  a  three-dimensional 
view  of  a  plane  wave  prior  to  and  after  encountering  atmospheric  turbulence.  Notice  that  one  of 
the  types  of  phase  distortion  introduced  by  the  atmosphere  is  an  overall  average  slope  or  “tilt"  to 
the  wave  front. 

For  the  simulation],  the  telescopic  imagingsystem  will  be  modeled  by  the  simple  lens  equivalent 
model  shown  in  Figure  2.4.  The  object  of  interest  is  a  star,  which  is  modeled  by  a  point  source 
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Figure  2.2.  Greyscale  diagram  for  wave  front  encountering  turbulence. 


at  infinite  distance  away.  The  atmospheric  turbulence  is  accounted  for  by  a  random  phase  screen, 
and  the  telescope  by  a  limiting  aperture  and  imaging  lens.  The  star  produces  a  planar  wave  front 
incident  on  the  phase  screen.  As  discussed  previously,  the  phase  screen  introduces  different  phase 
delays  along  the  wave  front  in  the  same  fashion  as  would  random  turbulence  to  a  planar  wave  front 
from  a  star.  This  degraded  wavefront  passes  through  the  aperture,  and  is  focused  by  the  lens  onto 
the  image  plane.  For  the  SRSH  hardware  shown  in  Figure  2.1,  the  wave  front  sensor  is  effectively 
located  in  the  aperture  plane  where  the  wave  front  is  sampled  prior  to  imaging. 

2.4  System  Transform  Function  (STF). 

2.4.1  Spatial  domain  image  formation  -  instantaneous  point  spread  function  (IPSF).  The 
goal  is  to  find  an  expression  for  the  impulse  response  of  the  optical  system  shown  in  Figure  2.4. 
Using  Fourier  system  concepts,  the  short  exposure  image  can  be  related  to  the  original  object  by  a 


2-3 


PLANE  WAVE  INCIDENT 
ON  EARTH'S  ATMOSPHERE 


Figure  2.3.  Effect  on  plane  wave  front  encountering  turbulence, 
convolution  relationship: 

i{x)  =  s(^*o(x),  (2.1) 

where  x  is  a  two-dimensional  vector  in  the  image  plane,  i  is  the  short  exposure  image,  o  is  the 
object  (star),  *  stands  for  a  convolution  operation,  and  s  is  a  function  called  the  instantaneous 
point  spread  function  (IPSF)  resulting  from  the  combined  effect  of  the  optics  and  the  atmosphere. 
For  a  single  star  in  the  field-of-view,  the  IPSF  represents  the  short  exposure  image. 
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Figure  2.4.  Lens  equivalent  model  for  telescopic  imaging  system. 

2.4.2  Frequency  domain  image  formation  -  system  transfer  function  (STF).  The  convo¬ 
lution  relationship  in  equation  (2.1)  can  be  Fourier  transformed  into  the  spatial  frequency  domain 
to  give  [11:4527]: 

/(^  =  5(p)0(^  (2.2) 

where  /,  S,  and  O  are  the  Fourier  transforms  of  «,  s,  and  o  respectively,  and  p  is  the  spatial 
frequency  vector.  For  the  case  of  a  star,  o{x)  is  effectively  a  point  source  an  infinite  distance,  and 
0{^  =  1.  This  reduces  equation  (2.2)  to: 

/(^  =  5(p1.  (2.3) 

The  function  S  in  equation  (2.2),  which  is  knqy^n  as  the  optical  transfer  function  (OTF),  is  deter¬ 
mined  by  the  expression  [11-4527],  [13:113-125]: 

_  (^(■'^^exp{j2;r<)>(AF^}]©[/’*(AF^exp{-;2n-<;i(AF’^}] 

■  [P(AFp)©P-(AFp)]  ’ 

where  P  is  the  pupil  funtion  which  defines  the  spatial  extent  of  the  aperture  of  the  telescope,  d) 
is  the  wave  front  phase  induced  by  the  phase  screen,  ♦  indicates  the  complex  conjugate.  F  is  the 
distance  from  the  pupil  plane  (aperture  in  Figure  2.4)  to  the  image  plane,  and  ®  indicates  an 
autocorrelation  operation. 
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Till'  owrall  ohjiTtivo  of  any  n'ronst ruction  tcrhniqnt*  is  to  go  from  tin'  turbnlcnci-  tlcgradod 
image  spectrmn.  /.  in  e(iualion  (2,2).  to  llie  host  estimate  achievable  (O)  of  the  original  object 
spectrum,  O.  An  I'stimale  of  the  original  object  is  then  obtained  by  an  invers<>  fourier  transform 
of  this  estimated  object  spectrum.  'I'he  approach  us»'d  for  the  SKSII  techiiiipie  involvi-s  u.sing  an 
estimate  of  the  OTF  to  invert  etpiation  (2.2).  The  problem  at  hand  is  to  in  some  w.ay  estimate 
Since  the  atmospheric  turbulence  effects  are  stvn  in  the  wave  front  pha.se  of  the  incident  wave,  wave 
front  si'iisor  mea.snrements  are  used  to  estimate  the  pha.se  in  tin'  pupil  so  that  the  estimate  for  the 
OTF  is  given  by  (1  l:4.'j27]: 


<■  [/’•(Arp)exp{-j2).0tAr^}] 

[P(AFe)eP*(AF^] 


where  .S  is  the  estimate  for  the  OTF,  and  <i  is  the  estimate  for  the  phase  screen  as  determined  from 
(he  wave  front  measiir  ineiits. 

Eqtiation  (2.2)  implies  the  original  object  spectrum  ran  be  recovered  by  simply  dividing 
through  by  the  OTF  estimate.  .S^,  The  problem  is  that  S,  as  an  estimate  of  5,  generally  will  not 
contain  all  the  information  in  >'  will  most  likely  contain  “zeros”  at  some  spatial  frequencies, 
indicating  lost  imaging  information.  A  much  better  approach  is  to  invert  equation  (2. .3)  through  a 
Wiener  filter  type  approach  so  as  to  minimize  the  mean  .square  error  of  the  estimate.  This  estimate 
of  the  object  spectrum  is  given  as  [I '.;d.V27]:  _ _ 


0(p)  =  lip) 


Wiener  type  filter 


(2.6) 


where  O  is  the  estimate  of  the  object  spectrum  and  t  is  a  parametric  constant  that  prevents 
sigularities  in  the  estimator.  Since  the  single  frame  OTF  estimate,  .S',  in  equation  (2.6)  will,  on 
average,  contain  zeros  at  some  spatial  frequencies,  the  estimate  of  the  object  spectrum  given  by- 
equation  (2.6)  will  be  incomplete  [M:2f)]. 


The  SRSII  approach  to  image  reconstruction  is  to  average  a  Wiener  type  filter  estimator 
over  many  frames  so  that  the  estimate  of  the  object  spectrum  (O)  includes  all  spatial  frequencies 
out  to  the  diffraction  limit  of  the  telescope.  There  are  two  possible  methods  of  averaging  for  the 
SRSH  technique,  and  both  will  be  considered  in  this  thesis.  For  the  first  method  considered,  the 
system  tran.sfer  function  (STf')  that  characterizes  the  SRSH  reconstruction  approach,  is  derived 
from  equation  (2.6)  by  averaging  the  numerator  and  denominator  to  give:  [15:2]: 

STFUf^  =  Q(p1  =  (2.7) 

where  STF\((>)  represents  an  ensemble  average  STF,  and  the  image  spectrum  /(p)  has  been  replaced 
by  S((i)  since,  as  previously  stated,  the  object  is  a  star.  The  notation  {•)  represents  averaging  over 
an  ensemble  of  a  large  number  of  short  exposure  realizath.ns  of  S(^  and  S(p).  This  equation  is 
used  to  characterize  the  performance  of  the  SRSH  technique  under  varying  seeing  conditions  for 
spatial  frequencies  out  to  the  diffraction  limit  of  the  telescope.  Alternately,  the  STF  could  be 
defined  by  averaging  the  ratio  of  numerator  and  denominator  to  give: 

=  (2.») 

where  STF2{^  represents  an  ensemble  average  STF. 

These  two  STF  estimators,  STFi  and  STF^,  represent  different  approaches  to  image  re¬ 
construction.  The  magnitude  of  a  STF  represents  the  spatial  frequency  response  of  the  SRSH 
technique  for  the  case  of  a  point  source  object.  The  results  for  both  SRSH  methods  are  presented 
in  Chapter  IV. 
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2.5  Assumptions  and  Defining  Equations  for  the  Thesis. 


This  thesis  is  a  follow-oii  to  earlier  work  tliat  involved  characterizing  the  ensemble  average 
performance  of  the  STFi  algorithm  for  varying  telescope  parameters  and  “seeing”  conditions  [lo]. 
Ihis  thesis  uses  limited  numbers  of  random  phase  screens  and  wave  front  measurements  to  estimate 
the  performance  of  STF\  and  STF^  algorithms.  Additionally,  this  thesis  presents  results  for  STFi 
and  ST Fj  signal-to-noise  ratios  (SNRs)  that  characterize  the  single  frame  STF  imaging  performance 
as  a  function  of  spatial  frequency  of  the  two  STFs.  The  assumptions  and  defining  equations  for  the 
simulation  results  are  presented  in  this  section. 

2.5.1  Assumptions. 

2. 5. 1.1  Required  sampling.  Simulation  runs  of  2000  frames  were  assumed  sufficient 
to  characterize  the  performance  of  the  estimator.  In  general,  the  SRSH  technique  requires  less 
frames  for  image  reconstruction. 

9 

2.5. 1.2  Assumed  turbulence  characteristics.  In  modeling  the  atmospheric  turbu¬ 
lence  by  a  phase  screen,  the  assumption  is  that  the  amplitude  of  the  incident  wave  front  is  not 
significantly  perturbed  by  the  turbulence.  This  turbulence  condition  is  commonly  referred  to  as 
near-field  turbulence.  For  astronomical  applications,  this  is  a  reasonable  approximation.  The  soft¬ 
ware  generated  random  screens  are  strictly  phase  screens  and  as  such  do  not  take  into  account  any 
amplitude  fluctuations. 

2.5.2  Defining  equations. 

2.5.2. 1  Phase  reconstruction.  In  a  practical  system,  slope  measurements  are  used  to 
reconstruct  an  estimate  of  the  incident  wave  front  phase  (!>{x)  through  a  linear  combination  of  zonal 
interpolation  functions.  For  this  thesis,  the  zonal  interpolation  functions  used  were  two-dimensional 
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gaussian  functions  specified  by: 


gi{x)  =  exp 


{ 


(2.9) 


where  Qi  is  the  *th  function  centered  at  x,-,  and  La,  is  the  radius  of  the  j  width  of  j,-.  Recall  from 
Figure  1.3  that  the  wave  front  sensor  divides  the  aperture  up  into  smaller  subapertures  over  which 
slope  measurements  are  taken.  When  reconstructing  the  wave  front  from  slope  measurements,  the 
zonal  interpolation  functions  are  located  at  the  corners  of  the  sensors  and  are  weighted  appropriately 
for  a  best  fit  of  the  slope  data. 

2. 5. 2. 2  Performance  related  equations  and  parameters.  Both  STFi  and  STF-j  and 
their  respective  variances  were  calculated  for  simulation  runs  consisting  of  2000  frames  of  simul¬ 
taneous  wave  front  sensor  measurements  and  short  exposure  images.  The  STF  results  are  plotted 
with  respect  to  a  normalized  spatial  frequency  defined  by; 


(2.10) 


where  ^  is  the  normalized  spatial  frequency,  p  is  the  unnormalized  spatial  frequency,  and  pjjjf  is 
the  magnitude  of  the  maximum  or  diffraction  limited  spatial  frequency. 

An  additional  way  to  characterize  the  STF  is  through  the  use  of  a  signal-to-noise  ratio  (SNR) 
measure.  A  single  frame  SNR  is  calculated  for  STF\  and  STFn  in  slightly  different  fashions.  For 
STFi,  we  note  that  the  purpose  of  the  denominator  is  ju.st  to  boost  the  high  spatial  frequencies. 
This  could  be  accomplished  by  a  deterministic  denominator.  All  the  Fourier  phase  correction  is  in 
the  numerator,  and  so  the  SNR  ratio  for  STFi  can  be  defined  by: 


SNRii^ 


\A’lium(^ 


(2.11) 
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where  SNRi  is  the  single  frame  SNR  for  STFi,  <’’num  variance  in  the  numerator  of  STFi, 

and  s//””"’  is  the  ith  realization  of  the  numerator  of  STFi.  For  STFo,  the  SNR  ratio  is  defined 
by: 

SNKM  =  .  (2.12) 

where  Sh'R2  is  the  single  frame  SNR  for  STFo  and  is  the  variance  of  STF2. 

The  SNRs  and  STFs  described  above  are  two  dimensional  equations.  For  this  thesis,  these 
two  dimensional  functions  are  reduced  to  one  dimensional  functions  for  presentation  purposes. 
The  reduction  is  accomplished  by  angularly  averaging  as  shown  in  the  following  reduction  of  the 
two-dimensional  function  g(pn)  to  the  single  dimensional  function 


9'{\on\)=^f\(\Pnie)dO, 


(2.13) 


where  |/51i|  and  0  are  the  polar  coordinate  repre.sentation  of  p\\.  These  angularly  averaged  results 
are  presented  in  Chapter  IV. 
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III.  SRSH  Computer  Simulation. 


3.1  Introduction. 

The  simulation  approach  to  generating  SRSH  transfer  functions  STF\  and  STF2  and  their 
respective  SNRs  avoids  the  rather  intensive  computations  required  for  an  analytical  approach.  In 
a  prior  thesis,  the  analytic  derivation  of  the  transfer  function  STFx  involved  a  four  dimensional 
numerical  integration  [16:3.6].  For  the  variance  calculations  required  in  the  SNR  analysis,  an  eight 
dimensional  integration  would  have  to  be  num<‘rically  evaluated. 

The  SRSH  simulation  is  a  modified  version  of  an  existing  adaptive  optics  FORTRAN  simu¬ 
lation  written  by  Dr.  Mik^  Roggemann  while  at  the  Phillips  Laboratory,  Kirtland  AFB,  NM.  The 
SRSH  simulation  uses  a  least  squares  phase  reconstruction  algorithm  to  produce  the  phase  front 
estimate  Additionally,  the  SRSH  simulation  u.ses  existing  tilt  and  piston  correction,  phase  .screen 

creation,  and  OTF  calculaUon  algorithms.  This  chapter  first  gives  a  description  of  simulation  steps 

[ 

taken  to  model  a  SRSH  e^qiped  telescope,  and  second  introduces  pertinent  analytical  equations 

used  in  the  simulation.  1 

I 

1 

3.2  Simulation  of  an  SRSH  Equiped  Telescope. 

The  SRSH  technique  involves  manipulation  of  the  telescope  OTF  (S)  and  its  estimate  (S) 

I 

over  many  frames.  The  sequence  of  operations  performed  by  the  SRSH  simulation  to  calculate 
OTFs  is  shown  in  Figure  3.1,  and  will  be  described  in  this  section.  The  simulation  applies  three 
different  templates  to  the  incident  wave  front  <l>.  A  WFS  template  and  an  associated  artificial  pupil 
P'(x)  are  used  for  WFS  meaisurements.  This  artificial  pupil  (P')  allows  for  WFS  subapertures 
with  corners  extending  outside  the  telescope  pupil  P(x)  to  still  be  used  for  WFS  measurements. 
The  third  template  is  the  telescope  pupil  P(x)  which  defines  the  aperture  of  the  telescope  and  is 
used  for  calculating  5  and  S.  These  templates  will  be  described  in  this  section  for  the  sequence  of 
simulation  steps  required  for  a  single  short  exposure  frame  calculation  of  S  and  S. 
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For  a  single  short  exposure  frame,  the  simulation  generates  a  random  phase  screen  or  wave 
front,  4>(x),  incident  on  the  telescope.  The  next  step  is  to  prepare  for  wave  front  spnsor(VVFS) 
slope  measurements  by  applying  an  artificial  pupil  P'  to  ^(F)  prior  to  application  of  the  WFS 
sensor  template.  This  is  done  because  the  slope  in  a  WFS  subaperture  is  actually  determined  by 
sampling  wave  front  at  points  along  the  WFS  subaperture  edges.  As  shown  in  Figure  3.2,  the 
WFS  has  subaperture  corners  which  extend  outside  the  telescope  pupil  P(x).  The  artificial  pupil 
defines  the  extent  of  4){x)  for  the  WFS  according  to  the  equation: 


where  D  is  the  telescope  or  aperture  diameter.  Both  pupils  are  shown  in  the  contour  plot  in 
Figure  3.3.  This  corrected  wave  front  is  now  used  to  determine  both  the  OTF  and  it’s  estimate  by 
the  two  separate  paths  shown  in  Figure  3.1. 

For  the  OTF  estimate  (5),  the  first  step  is  to  apply  the  wave  front  sensor  template  to  the 
corrected  wave  front.  The  result,  as  shown  in  Figure  3.1,  is  a  wave  front  sampled  along  the 


3-2 


defining  edges  of  the  wave  front  sensor  subaperlures.  The  reason  for  this  is  the  simulation  uses 
phase  differences  across  a  subaperture  to  calculate  the  slope  within  that  subaperture.  The  wave 
front  sensor  template  is  shown  in  Figure  3.2,  along  with  the  centers  of  the  Gaussian  interpolation 
functions.  The  estimate  of  the  phase  front  <l>  is  obtained  by  summing  the  spatial  respon.se  of  the 
interpolation  functions,  each  of  which  is  weighted  by  a  factor  determined  from  a  linear  combination 
of  the  slope  measurements.  The  telescope  puj)!!,  P,  is  applied  to  the  phase  estimate,  and  tilt  and 
piston  are  removed  from  the  phase  estimate  before  calculating  5. 

For  5(p),  the  first  step  is  to  apply  the  pupil  P  to  the  wave  front  previously  corrected  for 
WFS  measurements.  Next,  piston  and  tilt  (if  tilt  has  not  been  removed  prior  to  the  VVFS)  are 
removed.  The  OTF  is  then  calculated  from  the  resulting  corrected  phase  screen.  5(^  is  actually 
the  frequency  domain  result  of  a  single  frame  short  exposure  image  for  a  single  star.  Figure  3.4 
shows  the  simulated  image  plane  result.  The  result  is  similar  to  the  speckled  short  exposure  image 
shown  in  Chapter  I. 

3.3  Analytical  Equations  Used  in  the  Simulation. 

3.3.1  Flowchart  for  the  SRSH  simulation.  The  SRSH  computer  simulation  follows  the 
flowchart  shown  in  Figure  3.5.  The  three  main  sectipns  are  Setup,  Iterations,  and  Output.  A  brief 
overview  of  the  simulation  is  presented  in  this  section. 

3.3.2  Setup.  The  first  step  in  the  simulation  involves  specification  of  arrays  and  variables, 
and  initialization  of  simulation  parameters.  Data  is  read  in  to  specify  the  number  frames  .M  to 
process,  number  of  subapertures  accross  the  pupil  P,  telescope  diameter  D,  Fried  parameter  r„, 
flag  for  tilt  removal  (ITILT  in  Figure  3.5),  average  photon  count  per  subaperture  per  exposure 
N ,  a  small  parametric  constant  to  prevent  singularities  in  the  STF  estimator  e,  and  filenames  for 
simulation  output. 


CORMCTBD  WAVS  mOWT 
foa  WAVS  ntottr  8SMK>a 
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Figure  3.3.  Contour  plot  showing  size  of  actual  pupil,  P,  and  artificial  pupil,  P'. 

telescope  pupil  P  for  the  208  wave  front  sensor  subapertures  and  for  the  centers  of  the  241  gaussian 
reconstruction  functions  are  shown  in  Figure  3.2.  Gaussian  reconstruction  functions  are  specified 
by; 

jj(x)  =  exp| - (3.3) 

where  y,-  is  the  ith  function  centered  at  Xj,  x  is  the  position  vector  in  the  pupil,  and  Lat  is  the 
distance  between  adjacent  reconstruction  functions. 

In  addition  to  the  WFS  used  template  to  model  the  spatial  extent  of  the  subapertures,  the 
simulation  also  adds  in  random  noise  to  phcise  measurements.  The  accuracy  of  phase  reconstruction 
is  affected  by  both  the  signal  level  in  each  subaperture  of  the  wave  front  sensor,  and  the  strength 
of  atmospheric  turbulence.  Low  light  levels  lead  to  noisier  measurements.  Additionally,  the  Fried 
parameter,  r,,,  is  a  measure  of  turbulence  that  relates  to  the  spot  siae  and  so  must  also  be  considered. 
By  modeling  the  focused  spot  of  a  subaperture  lens  as  gaussian  in  shape,  the  standard  deviation 
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Figure  3.4.  Examp  le  of  simulation  short  exposure  image, 
for  phase  tilt  measurements  in  the  x  ot  y  directions  is  given  as  [18:1919]: 


/ 


where  <t„  (rad/m)  is  the  standard  deviation  of  the  tilt  measurement  and  is  a  dimensionless 
parameter  accounting  for  imperfections  in  the  subaperture  focal  plane  detector.  This  standard 
deviation  is  used  to  generate  random  numbers  which  are  added  into  phase  front  measurements  as 
noise. 

The  final  step  prior  to  entering  the  iterative  loop  of  the  simulation  is  to  derive  a  least- 
squares  recoiistructic  1  matrix,  R}‘ ,  which  is  used  for  converting  wave  front  slope  measurements  into 
weight  ngs  for  the  phase  interpolation  functions.  The  resulting  weightings  are  used  to  construct  a 
phase  estimate  by  summing  the  responses  of  the  241  appropriately  scaled  Gaussian  interpolation 
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functions  according  to  the  equation; 

241 

=  ^<'kah{^)  (3-5) 

t=i 

where  g*  is  the  Jbth  numbered  Gaussian  interpolation  function  located  at  a  corner  of  a  subaperture, 
and  cic  is  its  associated  weighting.  The  least  squares  approach  is  used  to  solve  for  the  required 
weighting  factors  by  inverting  the  linear  equation; 

S  =  HC,  (3.6) 

where  C  is  a  vector  containing  the  required  weightings  (ct),  5  is  a  vector  containing  slope  mea- 
suremen  (sf  and  //  is  a  Jacobian  matrix  whose  elements  are  defined  by  [19:457-458] 

[20:19-20],  For  the  particular  geometry  considered  here,  241  weightings  are  determined  from  416 
slope  measurements.  The  linear  system  described  by  equation  (3.6)  is  what  Tyson  describes  as 
“overdetermined  [4:237],’^  having  fewer  unknowns  (required  weightings)  than  linear  equations  to 
solve.  A  “best  fit”  to  the  data  (slope  measurements)  that  minimizes  the  mean  square  error  between 
the  left  and  right  sides  of  equation  (3.6)  gives  the  required  weightings  as: 

C  =  S,  (3.7) 

where  is  defined  eis  the  least  squares  reconstruction  matrix  used  to  map  slope  measurements 
into  actuator  commands  for  phase  estimate  construction. 

3.3.3  Iterations.  The  simulation  produces  results  for  M  short  exposure  frames.  Each  pass 
of  the  iteration  loop  involves  phase  screen  generation,  tilt  and  piston  correction,  phase  estimation, 
5  and  5  calculations,  and  summation  of  STF  results  with  those  from  previous  frames. 
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.  Phase  .irrrrn  grtirrafion.  Th<’  KL-FI'T  l<‘fhnique  is  used  to  genorato  random 

phaso  srrmis.  Tlio  tofliMi<iu<:‘  involves  the  creation  of  an  FF'l  screen  by  the  inverse  transform  of 
appropriately  filtered  complex  (Jaussian  random  variables.  'I'his  FFT  sceen  is  corrected  by  two 


low  order  KL  fimctions,  producing  the  complex  KL  FFT  screen  which  is  scaled  for  the 

turbulence  strength  (multiply  by  |p|  ).  'I’he  random  phase  screen  <i>  is  taken  from  the  real  or 
imaginary  part  of  ♦y jt-ki-  The  real  part  is  used  on  odd  and  the  imaginary  part  on  even  numbered 


frames.  The  details  of  this  technique  are  discussed  in  Appendix  A. 


3.S.3.2  Tilt  and  ptstoti  correction.  Next,  the  piston  error  of  the  random  phase 
screen,  0,  is  measured  in  the  telescope  pupil,  P,  and  removed  over  the  larger  pupil,  P',  as  follows 
(15:3-4]: 

-  J d'^xP[x)m  P'(x),  (3.8) 

where  d>'  is  the  piston  corrected  wave  front.  For  no  -tilt  simulations,  the  tilt  in  j;  and  y  directions 
is  calculated  over  the  pupil  P  as  follows  (21:253-255): 


J(Px{d„x}<i>{x)P[x) 
,fd^x{d„ .  *}2P(F) 


(3.9) 


where  („  is  the  tilt,  d„  is  a  unit  vector  in  the  x  or  y  direction.  In  preparation  for  wave  front  sensor 
measurements,  the  tilt  is  removed  over  P': 


=  (3.10) 

where  is  the  tilt  corrected  wave  front,  1,  and  ty  are  tilts,  and  i  and  J  are  unit  vectors  in  the  x 
and  y  directions. 

3. 3. 3. 3  Phase  reconstruction.  Since  most  of  the  results  for  this  thesis  are  for  the 
case  of  tilt  correction  prior  to  VVFS  measurements,  the  remaining  simulation  description  will  be 
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for  the  tilt  corrected  wave  front,  <t>"-  Slope  measurements  are  now  made  on  the  wave  front,  ip",  by 
applying  the  wave  front  sensor  template.  For  the  nth  siibaperture,  the  measured  slope  is  given  by 
[18:1914]: 

s„  =  J  J-x\V„(x)[^<l>"(x)  d„]  +  c,r„  (3.11) 

where  •  indicates  a  dot  product,  s„  is  the  slope  measurement  (radians  per  meter),  V  is  the  two- 
dimensional  spatial  gradient,  d„  is  a  unit  vector  in  the  direction  of  the  sensitivity  (i  or  y)  of  the 
nth  sensor,  and  a„  is  the  .slope  measurement  error  in  radians  per  meter.  Additionally,  IF„  (units 
of  m"*)  is  defined  as  a  weighting  function  for  the  nth  subaperture  such  that: 

tr„(if)  =  J  d^nv„{x)  =  1.  (3.12) 

The  slope  is  calculated  in  the  simulation  from  a  slightly  modified  version  of  equation  (3.11)  that 
expresses  the  integral  in  in  terms  of  phase  differences  along  the  edges  of  a  particular  subaperture. 

To  get  the  phase  difference  version  for  slope,  equation  (3.11)  is  first  rewritten  using  integration 
by  parts  to  give  [18:1914]: 

8„  =  -  J  d^x<p"{x)[VW„{x)  ■  d„]  +  a„,  (3.13) 


where  the  integral  is  now  expressed  in  terms  of  the  gradient  of  the  weighting  function  for  the  nth 
subaperture,  VWn(-^-  For  the  subaperture  shown  in  Figure  3.6,  this  gradient  can  be  expressed  as: 


VVF„(x)  = 


u.sed  to  calculate  slope  in  y  direction 

. . . . .  ^ 

/  dx{\V„{x,yt)6(x,y-yi)  -  W„{x,yr)6{x,y  -  yr)}  + 

Jx\ 
yyi 

/  dy{W„(xr,y)(i(x  -  Xr,y)  -  VF„(x,,y)(5(x  -  j,,y)}, 

*'yr 


used  to  calculate  slope  in  x  direction 


L  =  Subaperture 
side.  5  pixels 


(x,.y,l 

•  • 


•  • 

•  • 

'*r-V 

Figure  3.6.  Location  of  points  on  nth  subaperture. 

where  6(x,y)  is  an  impulse  function,  and  (xi,yi)  is  the  coordinate  of  the  upper  left  corner  and 
(^riJ/r)  the  lower  right  corner  of  the  nth  subaperture  as  shown  in  Figure  3.6.  Using  the  “sifting” 
property  of  an  impulse  function  [22:56-57]  equation  (3.13)  can  be  rewritten  in  terms  of  line  integral 
phase  differences  along  the  subaperture  edge.  For  example,  the  y  slope  for  the  nth  subaperture  is 
given  by: 

»n  =  /  d^<l>"{x,yi)W„(x,yi)-  f  dx (j>(x,yr)W„{x,yr)  +  a’ft,  (3.15) 

Jtt  Jxi 

where  sjj  is  the  y  slope  for  the  nth  subaperture,  and  ajj  is  a  the  additive  random  noise  term  of  zero 
mean  and  standard  deviation  determined  by  Equation  (3.4). 

The  reconstruction  matrix,  R^’ ,  is  used  to  map  the  vector  of  slope  measurements  onto  weight¬ 
ing  commands  for  phase  reconstruction.  The  vector  of  weighting  commands  is  determined  by: 

C  =  R'’S,  (3.16) 


where  C  is  the  vector  containing  241  actuator  commands  for  scaling  241  gaussian  functions,  S 
is  the  vector  containing  x  and  y  slope  measurements  (416,  two  measurements  for  each  of  208 
subapertures),  and  Rf'  is  the  reconstruction  matrix  (241  rows  by  2  x  208  =  416  columns).  The 
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phase  est  imate  is  then  constructed  by  summing  the  scaled  gaussian  function  responses  at  each  point 
in  the  pupil  to  give: 


<^{x)  = 


(3.17) 


where  0  is  the  phase  estimate,  c*  is  the  itth  weighting  command  or  scaling  for  the  i-th  gaussian 
function  located  at  £,■. 

3. 3. 3. 4  Calculating  centered  OTFs  for  S  and  S.  For  the  SRSH  technique,  the 
simulation  calculates  centered  OTFs  for  5  and  5.  Centering  the  OTF  is  accomplished  by  removing 
tilt  from  the  phase  fronts.  Since  this  discussion  is  for  tilt  removed  prior  to  WFS  measurements, 
<t>"  is  the  corrected  wave  front  for  calculating  5,  and  as  such  already  has  tilt  removed.  For  5,  the 
phase  estimate  0  does  require  tilt  correction.  Additionally, the  telescope  pupil  P  is  applied  to  <*>" 
and  ^  in  preparation  for  OTF  calculations. 

5(^  is  calculated  directly  from  a  forward  Fourier  transform  of  the  magnitude  squared  of  the 
coherent  impulse  response  [13:114]: 


(3.18) 


where  P'{-}  indicates  a  forward  Fourier  transform,  and  A,  is  the  coherent  impulse  response  for  the 
Jth  frame  as  defined  by  [13:111]: 


hi{x)  =  :r{P(AF^exp{-2:r<(AF^)}  ,  (3.19) 


where  A  is  the  wavelength,  and  F  is  the  telescope  focal  length.  The  estimate  of  the  OTF,  5,-,  is 
obtained  if  is  used  instead  of  0"  in  equation  (3.19). 


Next,  STF  quantities  are  calculated.  For  STF\,  the  single  frame  numerator  and  denominator 
are  calculated  by: 

(3.20) 


and; 


s</^'^"(^  =  |S.(p)p  +  c, 


(3.21) 


where  <  is  an  arbitrarily  small  number  to  prevent  divide  by  zero,  st  is  the  ith  frame  numerator, 
and  is  the  ith  frame  denominator.  For  STF^,  the  single  frame  result  is  calculated  as: 


(3.22) 


where  stfi  is  ith  frame  result. 


3.3. 4  Results  of  simulation  run.  The  STF’s  are  calculated  by  averaging  the  single  frames 
results  over  M  frames.  So,  STFI  is  calculated  as  follows: 


Finally,  STF2  is  found  by; 


STF2(pi  =  {s</.  (p)) , 


(3.23) 


(3.24) 


where  stfi(^  is  the  ith  single  frame  result. 

The  SNR  for  a  single  frame  of  the  image  can  be  defined  in  the  spatial  frequency  domain  as 


[19:452]: 


SNR(p)  = 


(3.25) 
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0 


where  7  is  the  Fourier  transform  of  the  image  and  tr]{p)  is  the  image  variance.  For  SNR  calculations, 
the  variance  is  calculated  for  functions  according  to  the  equation  [23:244]: 


(3.26) 


where  <rj(^  is  the  variance  of  the  random  process  /(^,  and  /i(^  is  j'th  frame  realization  of  the 
function.  Instead  of  an  image,  this  thesis  is  interested  in  the  SNR  of  the  STF  estimators.  For 
STFi,  the  SNR  is  defined  in  terms  of  the  comple.x  numerator,  sifP^^.  The  single  frame  SNR  for 


STFi  is  defined  by: 


SNR,(^ 


\/‘^num(^ 


(3.27) 


where  <rnum(^  ‘s  the  variance  of  the  STFi(p)  numerator.  For  STF^,  the  single  frame  SNR  is 
defined  in  a  more  straightforward  fashion  as: 


(3.28) 


where  is  the  variance  of  The  results  for  SNRs  and  STFs  are  presented  in  the  next 


IV.  Simulation  Result.^. 


4.1  Introduction. 

This  chapter  presents  simulation  results  of  the  system  transfer  function  (STF)  and  associated 
signal-to-noise  ratio  (SNR)  measures  defined  in  Chapter  II.  This  chapter  will  begin  with  a  review 
of  previous  analytic  results  for  the  SRSH  technique  found  in  Reference  [24].  These  analytic  results 
characterize  the  STF,  and  as  such,  serve  as  the  basis  for  this  thesis  effort.  General  analytic  STF 
trends  will  be  used  for  discussing  the  validity  of  simulation  STF  results. 

Next,  STFi  and  STF2  re.sults  are  presented.  The  simulation  results  were  generated  using 
2000  randomly  generated  phase  screens  and  the  conditions  summarized  in  Table  4.1.  The  telescope 
diameter  (D)  was  set  at  two  meters.  Also,  the  parametric  constant  (f)  to  prevent  divide  by  zero 
in  STF  calculations  was  set  at  1  x  10~®.  Only  selected  results  from  these  runs  are  presented  in 
this  thesis.  Individual  STF  simulation  runs  are  characterized  by  a  subaperture  light  level  (average 
photons  per  subaperture  per  short  exposure  image,  N),  possibility  of  tilt  correction,  and  ratio 
of  center-to-center  subaperture  separation  to  Fried  parameter,  This  ratio,  defines  the 
adequacy  of  the  wave  front  sampling  in  detecting  the  turbulence  characteristics  of  the  incident 
wave  front.  Although  correcting  tilt  prior  to  the  WFS  is  optional,  the  majority  of  the  results 
discussed  will  be  for  the  case  of  tilt  corrected  phase  fronts.  Not  removing  tilt  leads  to  a  higher 
error  in  the  phase  estimate.  The  effect  of  no-tilt  removal  will  be  briefly  examined  for  STFi  ■ 

This  chapter  also  presents  an  analysis  of  the  non-zero  mean  bias  terms  in  the  numerator  and 
denominator  of  the  STF.  The  results  indicate  the  OTF  estimate,  S,  contains  a  negative  bias.  The 
bias  problem  is  particularly  significant  under  the  combined  conditions  of  low  incident  light  levels 
and  good  WFS  sampling  (sampling  ratio  ^  less  than  one). 

The  final  simulation  results  presented  are  selected  frequency  domain  SNRs  which  characterize, 
as  a  function  of  spatial  frequency,  the  single  frame  STF  image  reconstruction  performance.  The 
SNRs  can  be  used  to  determine  the  number  of  short  exposure  frames  of  data  to  obtain  a  required 
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PARAMETER 


VALUES  USED 


PHOTONS  PER  SUDAPERTIIRE 

PER  SHORT  EXPOSliRE  IMAGE 

10.  20.  50.  100.  200.  400. 

800.  1000 

RATIO  OF  SUBAPERTURE  SIZE 

TO  FRIED  PARAMETER.  L/r^ 

0.25,  0.5,  I.O.  2.0 

PHASE  CORRECTIONS  PRIOR 

TO  WAVE  FRONT  SENSOR 

NONE.  OR  CORRECT  TILT 

CONTROL  LAW  USED  FOR 
CONSTRUCTING  ESTIMATE 

OFTHE  PHASE  i 

MINIMUM  MEAN 

SQUARE  ERROR 

Table  4.1.  Scenarios  considered  for  thesis  research. 

SNR  level  out  to  a  given  spatial  frequency.  The  analytic  con.sideration  of  Reference  [24]  did  not 
include  any  SNR  results  due  to  the  required  eight-dimensional  numeric  integration  for  the  STF 
variance.  The  ratio  is  presented  to  discuss  the  comparative  advantages  of  STF\  and  STF2. 

4.2  Summary  of  Previous  Analytic  Results  [24:15-16]. 

The  analytic  results  presented  in  this  section  are  for  a  numeric  integration  of  a  derived  four 
dimensional  integral  expression  for  the  STF  in  terms  of  the  OTF  (5)  and  its  estimate  (S)  [24:10]. 
Due  to  the  complexity  of  this  integration,  the  analytic  results  are  for  case  of  a  simple  square 
aperture  telescope  and  a  small  nur.iber  of  WFS  subapertures.  The  WFS  used  has  a  maximum  of 
five  subapertures  accross  the  telescope  pupil  as  shown  in  Figure  4.1.  Analytic  results  are  presented 
for  variations  in  both  the  average  photon  count  per  subaperture  per  short  exposure  image  N, 
and  the  WFS  sampling  ratio,  Analytic  result  trends  are  used  for  discussing  the  validity  of 
simulation  STF  results  in  the  next  section.  The  comparison  is  not  rigorous  since  analytic  STF 
results  are  for  a  different  phase  estimate  construction  control  law  (minimum  variance),  and  are  for 
a  square  telescope  pupil. 


The  analytic  STF  is  characterized  for  variations  in  the  both  the  the  average  photon  count 
per  subaperture  per  short  exposure  image  (A')  and  the  sampling  ratio  (^).  In  general,  WFS 
measurements  improve  with  increasing  signal  strength  of  the  incident  light  as  specified  by  increasing 
N.  The  other  important  parameter  concerns  how  dense  the  WFS  nrieasurements  arc  compared  to  the 
Fried  parameter  To.  The  slopes  are  measured  in  WFS  subapertures  with  center-to-center  spacings  of 
La,  •  The  ratio  of  ^  is  used  to  compare  the  slope  measurement  spacing  to  the  turbulence  conditions 
indicated  by  the  Fried  parameter.  Smaller  ^  values  signify  more  dense  WFS  measurements  with 
respect  to  the  size  of  turbulence  induced  distortions  on  the  incident  wave  front  since  r„  characterizes 
the  approximate  size  of  significant  wave  front  distortions.  If  the  sampling  ratio  is  less  than  unity, 
most  of  the  turbulence  characteristics  along  the  incident  wave  front  are  being  accurately  measured 
since  La,  is  smaller  than  Vg.  Poor  sampling  is  characterized  by  sampling  ratios  greater  than  unity. 
For  the  analytic  results,  N  ranges  from  10  to  400,  and  from  0.1  to  2. 

The  analytic  results,  which  are  from  a  paper  by  VonNiederhausen  and  Welsh,  are  summarized 
in  Figures  4.2  and  4.3  [24:15-16].  The  STF  results  are  plotted  versus  a  one-dimensional  normal¬ 
ized  spatial  frequency,  where  fi  represents  the  x-component  of  the  two-dimensional  spatial 

frequency  vector  p,  and  F  is  the  focal  length  of  the  telescope. 

The  STF  is  a  measure  of  the  spatial  frequency  response  of  the  SRSH  reconstruction  technique. 
Ideally,  the  STF  would  be  unity  out  to  the  diffraction  limit  of  the  telescope.  For  a  sampling  ratio 
of  unit)^,  Figure  4.2  shows  how  the  SRSH  technique  performance  improves  toward  ideal  behavior 
with  increasing  light  intensity  (N  increases).  The  curves  asymptotically  approach  an  upper  limit 
im'  'sediby  the  sampling  ratio. 

Figure  4.3  illustrates  the  positive  result  on  STF  performance  obtained  by  decreasing  the 
sampling  ratio.  For  these  results,  the  average  photon  count  is  held  at  N  =  50,  while  the  sampling 
ratio  ranges  from  0.1  to  2.  The  STF  curves  indicate  that  the  technique  is  extremely  sensitive  to 
the  sampling  ratio. 
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Figure  4.1.  Telescope  aperture  and  wave  front  sensor  geometry  for  analytic  resi  i*s  [24:13]. 

In  summary,  both  light  intensity  level  (.A^)  and  sampling  ratio  (^)  are  importaat  parameters 
for  characterizing  expected  STF  image  estimator  performance.  The  simulation  STF  ';ehavior  will 
now  be  examined  for  variations  in  these  two  parameters. 


Figure  4.2.  System  transfer  function  (STF)  results  for  square  aperture  and  WFS  shown  in  Fig¬ 
ure  4.1  and  ^  =  1  [24:15]. 
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Figure  4.3.  System  transfer  function  (STF)  results  for  square  aperture  and  WFS  shown  in  Fig¬ 
ure  4.1  and  iV  =  50  [24:16]. 

4-3  STF  Results. 

The  simulation  STF  results  will  be  examined  in  terms  of  the  general  trends  noted  for  the 
analytic  STF  results  in  the  previous  section  [16:4.8-4.15].  Specifically,  STF  curves  should  approach 
asymptotic  behavior  with  increasing  photon  count.  This  asymptotic  value  must  be  equal  to  or 
less  than  unity.  The  STF  curves  should  also  improve  with  an  improved  sampling  ratio,  again 
approaching  but  not  exceeding  unity.  Deviations  from  the  trends  established  by  the  analytic  results 
will  be  explained  in  the  next  section  which  investigates  a  negative  noise  bias  problem  originating 
with  the  OTF  estimate,  S.  This  section  presents  a  summary  of  STF  results  for  the  case  of  tilt 
removed  prior  to  WFS  measurements.  This  is  the  baseline  case  for  this  thesis.  To  examine  the 
benefit  of  tilt  removal,  ST Fi  results  for  tilt  removed  and  tilt  not  removed  are  presented  for  a  photon 
count  of  100,  and  sampling  ratios  varying  from  0.25  to  2. 

Figures  4.4  and  4.5  show  the  tilt  removed  STF  curves  for  photon  counts  of  10  and  800,  with 
sampling  ratios  varying  from  0.25  to  2.  STFi  and  STFj  are  characterized  by  non-zero  valued  image 
reconstruction  out  to  the  diffraction  limit  of  the  telescope  for  all  sampling  conditions  considered. 
Except  for  a  sampling  ratio  of  0.25,  both  ST Fj  and  ST F2  show  improved  performance  for  increasing 
the  photon  count  while  holding  the  sampling  ratio  fixed.  Improving  the  sampling  ratio  also  generally 
improves  the  STF. 
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The  problem  with  the  simulation  results  is  that  both  STFi  and  STF^  curves  exceed  unity 
at  some  point  for  most  photon  counts  and  all  sampling  ratios  considered.  This  deviation  from  the 
previously  predicted  results  [24]  is  due  to  a  non-zero  bias  in  the  OTF  estimate  5.  The  bias  appears 
to  get  worse  as  the  sampling  varies  from  poor  >  1.0)  to  good  <  1.0)  for  both  STF 
estimator  methods.  For  each  sampling  ratio,  the  problem  gets  successively  worst  for  decreasing 
photon  counts,  indicating  the  S  bias  is  increased  with  increased  noise  in  the  slope  measurement.^ 
since  the  standard  deviation  of  the  additive  noise  in  slope  measurements  is  proportional  to 
The  best  example  of  this  trend  is  shown  in  Figures  4.4  and  4.5  for  sl  sampling  ratio  of  ^  =  0.25, 
where  the  STF  curve  for  the  low  photon  count  of  10  stays  well  above  unity,  except  for  spatial 
frequencies  out  near  the  diffraction  limit  (p„  =>  1).  The  STF  curve  for  the  high  photon  count  of 
800  is  closer  to  unity  for  out  to  the  diffraction  limit.  These  results  imply  the  spatial  frequencies 
of  the  image  are  being  amplified  to  give  better  than  the  theoretical  maximum.  This  is  clearly  an 
unphysical  result  which  has  not  been  previously  noted  in  the  literature.  The  next  section  traces 
this  result  to  the  bias  in  5.  The  bias  is  shown  to  get  progressively  worse  as  sampling  improves 
(smaller  sampling  ratio)  and  the  photon  count  remains  low. 

Figure  4.6  shows  the  effect  tilt  removal  has  on  STFi  for  a  photon  count  of  100,  and  varying 
sampling  conditions.  At  unity  and  low  sampling  ratios,  removing  tilt  prior  to  the  WFS  boosts  the 
STFi  estimator  curves  slightly  in  magnitude.  This  effect  is  most  noticeable  at  spatial  frequencies 
near  the  diffraction  limit.  As  WFS  sampling  is  improved  (sampling  ratio  drops  below  one),  removing 
tilt  gives  some  improvement  at  higher  spatial  frequencies,  with  negligible  results  at  lower  spatial 
frequencies. 

Noise  characteristics  of  the  simulation  STF  estimator. 

The  goal  of  this  section  is  to  characterize  the  bias  terms  in  the  STF.  To  investigate  the 
problem,  the  single  frame  OTF  (5)  is  expressed  as  a  random  complex  quantity  and  the  OTF 
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Figure  4.4.  Tilt  removed  STF\  curves  for  different  wave  front  sampling  conditions  as  specified  by 

Ltx 

r.  ■ 


estimate  (5)  as  noisy  version  of  5.  The  noise,  M,  is  then  determined  from  the  difference  between 


S  and  S  (5  —  S  =  M).  The  complex  expressions  for  S  and  5  are  used  to  derive  expressions  for  the 


numerator  and  denominator  of  the  STF  estimator.  An  investigation  of  the  noise  terms  in  both  the 


numerator  and  denominator  shows  the  source  of  the  problem  is  the  negative  mean  of  the  real  part 
of  .V  in  S.  An  attempt  will  be  made  to  identify,  from  the  bias  results,  spatial  frequency  ranges 


which  yield  valid  simulation  STF  results  (small  bias  regions). 


The  first  goal  of  this  analysis  is  to  isolate  the  complex  valued  additive  noise  term  contained 


in  S.  To  do  this,  5  is  first  expressed  as  a  complex  quantity; 


(4.1) 


where  Sut  and  Sx„  are  the  real  and  imaginary  parts  of  the  5,  and  j  indicates  the  imaginary  part. 
An  expression  for  S  is  found  by  adding  real  and  imaginary  noise  compcuents  to  S  in  equation  (4.1) 


to  give: 


+  A/We(^  +  j  ,  (4.2) 


I  ■ 


i 
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Figure  4,5,  Tilt  removed  ST  A  curves  for  different  wave  front  sampling  conditions  as  specified  by 

Lai. 

* 

where  and  A'jm  are  the  real  and  imaginary  parts  of  the  additive  noise  in  5.  To  isolate  the 
noise  in  5,  equation  (4.1)  is  subtracted  from  equation  (4.2)  to  give: 

5(p-)  -  S(p)  =  +  iA^xm(p)  =  A'(p)  (4.3) 

where  N  is  the  complex  noise  portion  of  5.  The  bias  problem  originates  from  a  negative  mean  real 
noise  term  in  A : 

{Ar{^)  =  {ynei^)<0,  (4.4) 

with  the  imaginary  noise  part,  jVjm,  having  zero  mean. 

The  next  step  in  the  noise  analysis  is  to  derive  complex  expressions  for  the  numerator  (5S* ) 
and  denominator  ((5p)  terms  of  the  STF,  and  average  the  expressions  to  determine  significant 
noise  terms.  The  complex  numerator  is  expressed  as; 

5(^5‘(pl  =  +  +  + 
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Figiiro  STF\  rurve.s  for  different  wave  front  sampling  conditions  a,s  specified  by  and  for 

tilt  and  no  tilt  removed. 


(4.5) 


As  will  be  shown  later,  the  only  significant  noi.sc  term  contributing  to  an  STF  numerator  bias  is 
So,  (5(p)5*)  can  be  written  as: 


{S{^S-(^)  =  (S|.)  +  {SnAf^A^nAf^)  +  (5|„.).  (d.e) 


The  complex  STF  denominator  term  (I.S|^)  is  derived  by  taking  the  stpiared  magnitude  of  equa¬ 


tion  (4.2)  to  give: 


=  •SL(f7)  +  -b.v^,(p)  +  SlJfT)  +  2Si,Air).Vi,„{p)  +  .Vi„.(p).  (4.7) 


The  significant  hi, as  terms  in  the  denoiiiinator  are  2.‘^KfAK,  and  A'^,  so  that  {|.S(/5)|-)  ran  be 


written  as: 


(l*‘'’(f7)P)  =  +  nSnApiXnAPl)  +  +  (.‘'1,„(Ar)). 


(4,8) 


/ 


A  comparison  of  equat  ions  (4.6)  and  (4.8)  leads  to  a  convenient  figure  of  merit  R  ^n-biased  STF 
performance  as  those  spatial  frequencies  over  which  the  numerator  hi;ts  tf'rm  is  eq;  or  less  than 
the  denominator  bias  terms.  This  approach  yields  the  following  ineipiality; 

<‘2{SR,{pj.\'T,4ir])  +  (A^h(^)  (4-9) 

> - ^  — I - ' 

numerator  bias  term  denominator  bias  term 

The  range  for  a  non-biased  estimator  is  easily  determined  by  plotting 
and  noting  where  the  ratio  goes  above  0.5. 

The  mean  valued  results  (jV'r,)  and  {jVjm)  are  plotted  in  Figure  4.7  for  a  photon  count  of 
10,  and  a  sampling  ratio  of  0.25.  Notice  that  only  the  real  part  (Afite)  is  significant.  The  results 
show  the  real  noise  term  is  negatively  biased.  The  imaginary  part  is  many  orders  of  magnitude 
less,  appearing  as  a  straight  line  in  Figure  4.7. 

The  noi.se  bias  problem  in  the  denominator  is  shown  in  Figure  4.8.  Since  the  real  part  of 
the  complex  noise  is  negatively  biased,  the  2{Sr,A/k<;)  noise  term  in  equation  (4.7)  is  the  most 
significant  since  it  involves  the  cross  product  between  the  negatively  biased  noise  term  and  the 
larger  and  positive  real  part  of  S  (.SRe)-  2(Sr,A'r,)  is  easily  identified  as  the  noi.se  term  most 
responsible  for  decreasing  the  magnitude  of  the  STF’  denominator.  This  is  easily  seen  in  Figure  4.8, 
where  ‘2{:.'iieAfne)  is  plotted  versus  all  denomoinator  noise  terms  (terms  with  or  Aj,,,  in 
them)  in  equation  (4.7).  The  (vVr,)  term  is  the  dominant  term  that  causes  the  denominator  bias 
to  decrease  slightly. 

For  the  numerator,  the  complex  portion  of  equation  (4.5)  essentially  averages  to  zero.  Of  the 
real  numerator  noise  terms,  only  the  (5R,A/Rf)  term  is  significant  as  shown  in  Figure  4.9.  The 
curve  for  all  real  numerator  noise  terms  and  just  {.‘'r,Ar,)  are  almost  overlapping.  A  comparison 
of  Figures  4.9  and  4.8  .shows  the  numerator  noise  terms  contribute  slightly  over  half  the  negative 
bias  as  the  denominator  noise  terms  for  these  photon  and  sampling  conditions.  Sinci-  lioth  noise 
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Figure  4.7.  OTF  estimate  (S')  real  and  imaginary  noise  content,  at  a  samping  ratio  of  ^  =  0.25 
and  a  photon  count  of  10. 

term  contributions  are  negative,  this  difference  between  the  numerator  and  denominator  is  what 
makes  the  STF  estimate  go  above  one. 

The  last  part  of  the  noise  analysis  is  to  define  valid  unbiased  regions  for  simulation  STF 
results.  The  ratio  is  plotted  in  Figure  4.10.  Where  this  ratio  goes  above  0.5 

is  where  the  denominator  bias  is  less  than  the  numerator  bias.  This  will  define  the  non-biased 
estimator  for  this  thesis.  For  a  sampling  ratio  of  ^  =  2,  the  criterion  used  indicates  the  STF  will 
be  unbiased  from  a  normalized  spatial  frequenccy  of  0.05  out  to  the  diffraction  limit.  At  a  unity 
sampling  ratio  =  1),  the  unbiased  STF  is  defined  from  approximately  a  normalized  spatial 
frequency  of  0.09  out  to  the  diffraction  limit.  At  a  sampling  ratio  of  ^  =  0.5,  the  estimator  is 
unbiased  from  about  0.3  on  out  to  the  diffraction  limit.  This  sampling  ratio  is  right  on  the  edge 
of  being  biased  out  to  the  diffraction  limit  as  seen  by  how  clo.se  the  ratio  stays  to 

the  value  of  0.5.  For  the  sampling  ratio  of  0.25,  the  STF  essentially  remains  biased  out  to  near  the 
diffraction  limit. 

Thesis  conclusions  and  recommendations  are  discussed  in  the  next  chapter. 
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Figure  4.8.  I  STF  denominator  noise  terms  versus  just  2STieJ^Tie,  at  a  samping  ratio  of  ^  =  0.25 
and  a  photon  count  of  10. 

■ 

4-5  SNR  Results. 

The  second  goal  of  this  thesis  is  to  present  single  frame  SNR  results  for  the  two  STF’s  under 

i 

varying  seeing  conditions.  These  single  frame  SNR  estimates  characterize,  as  a  function  of  spatial 
frequency,  the  required  number  of  fr’’.mes  of  short  exposure  images  and  WFS  measurements  to 
overcome  iniage  reconstruction  noise.  Only  the  results  for  SN R\  for  the  case  of  tilt  removed  prior 
to  WFS  measurements  are  presented  since  the  same  trends  were  noted  whether  tilt  was  or  was 
not  removed.  SNR2  behavior  is  represented  by  curves  for  the  ratio  snHj  •  SNR  results  are 
presented  for  the  sampling  ratio  of  0.25  due  to  the  previously  discussed  noise  bias  problem  for  the 
STF  estimator  at  that  sample  ratio.  The  results  are  presented  in  subsections  for  the  sampling 
ratio  conditions  of  unity  sampling  =  1.0),  undersampling  =  2.0),  and  oversampling 
(^  =  0.5). 

4.5.1  Unity  sampling,  ^  =  1.0.  Estimated  single  frame  5,VRi  results  for  unity  sampling 
are  shown  in  Figure  4.11.  The  SNR  approaches  asymptotic  behavior  above  photon  counts  of  50, 
indicating  STFi  is  increasingly  susceptible  to  noise  at  or  below  this  photon  count.  From  the  ratio 


\ 
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Figure  4.9.  STF  numerator  real  noise  terms  versus  jW  S■K.t■^^ne,  at  a  sampling  ratio  of  ^  =  0.2.5 
and  a  photon  count  of  10.  | 

curves  shown  in  Figure  4.12,  it  is  evident  STFi  has  better  SNR  characteristics  than  STF2 
as  indicated  by  the  ratio  value  exceeding  one  for  nortnalized  spatial  frequencies  above  roughly  0.1. 
There  is  roughly  a  factor  of  2  difference  between  SNRi  and  SNR2  almost  out  to  the  diffraction 
limit  of  the  telescope  (normalized  spatial  frequency  of  1.0).  One  interesting  point  to  note  is  that 
the  advantage  of  STFi  over  STF2,  as  indicated  by  th^  SNR  ratio,  decreases  with  increasing  photon 
count.  This  result  indicates  that  STF2  is  more  susceptible  to  the  the  higher  WFS  noise  encountered 
for  lower  photon  counts.  | 


4.5.2  Under  sampling,  =  2.0.  Figure  4.13  shows  results  for  SNRi  for  the  case 
of  under  sampling.  The  undersampled  SNR  approaches  asymptotic  behavior  at  a  photon  count 
exceeding  appoximately  100.  The  asymptotic  behavior  for  the  undersampled  .SNRi  occurs  at 
higher  photon  counts  than  the  previously  discussed  unity  sampling  case.  The  undersampled  WFS 
requires  higher  signal  levels  to  overcome  the  noise  inherent  in  poorer  sampling.  The  ratio 
curves  in  Figure  4.14  show  that,  except  for  low  photon  counts  and  high  spatial  frequencies,  STFi 
is  the  better  method  for  image  reconstruction  than  STF2.  The  advantage  of  SNRi  over  SNR2 
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Figure  4.10.  Ratio  curves  for  fof  <*  photon  count  of  10,  and  varying  sampling. 

begins  at  around  a  normalized  spatial  frequency  of  0.05,  with  approaching  3  for  iihoton 

counts  above  10.  For  the  low  photon  count  of  10  and  normalized  .spatial  frequencies  approaching 
unity,  the  two  STF  SNRs  are  roughly  equal.  The  STFi  image  estimator  approach  is  increasingly 
preferred  over  ST as  the  photon  counts  increases. 

4.5.3  Over  sampling,  ^  =  0.5.  Only  SNR  results  for  a  sampling  ratio  of  0.5  are 
presented  due  to  the  significant  STF  noise  bias  .STF  at  a  sampling  ratio  of  0.25.  The  SNR  results 
for  a  sampling  ratio  of  0.5  are  shown  in  Figure  4.15  for  SNlii.  As  expected,  the  SNRi  curves 
are  higher  than  the  previously  considered  under  =  2)  and  unity  =  1)  sampling  cases. 
The  better  results  are  due  to  a  more  accurate  phase  estimate,  leading  to  better  STF  estimator 
performance,  particularly  at  higher  spatial  frequencies.  Additionally,  the  advantage  of  the  STFi 
over  the  STFj  estimator  decreases  with  increasing  photon  count,  as  evidenced  by  curves 

approaching  unity  for  higher  photon  counts  in  Figure  4.16.  What  this  essentially  means  is  that  the 
STFi  and  STF2  estimators  are  equally  good  for  image  reconstruction  for  an  adequate  signal,  and 
a  good  sampling  ratio. 
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V.  Conclusions  and  Recommendations. 

5.1  Conclusions. 

This  thesis  is  a  follow-on  to  a  previous  analytic  investigation  of  the  STF  for  the  SRSII  im¬ 
age  reconstruction  technique.  The  thesis  objective  was  to  characterize  two  different  estimator 
approaches,  ST Fi  and  ST Fq  ,  and  to  compare  the  signal-to-noise  ratios  (SNRs)  to  determine  which 
technique  was  most  advantageous  for  image  reconstruction.  The  simulation  SNRs  were  also  included 
to  provide  a  means  to  determine  frame  requirements  to  achieve  a  desired  image  reconstruction  SNR. 
The  mo.st  significant  result  for  this  thesis,  however,  was  the  isolation  of  a  noise  bias  in  the  OTF 
estimate,  S. 

Ignoring  for  a  moment  the  bias  problem  encounted  in  the  simulation,  the  non-biased  results 
followed  trends  established  by  earlier  analytic  results  [2d].  In  general,  the  STF  can  be  improved 
with  increased  signal  strength  (characterized  by  the  average  photon  count  per  subaperture  per  .short 
exposure  image, A^)  up  to  an  asymptotic  limit  determined  by  the  wave  front  sensor  (WFS)  sampling 
ratio  {^).  Some  additional  are  the  STF  is  extremely  sensitive  to  variations  in  the  sampling  ratio 
for  a  fixed  signal  strength,  and  can  be  significantly  boosted  by  improving  the  sampling  ratio.  Finally, 
tilt  removal  prior  to  WFS  measurements  made  a  slight  improvement  on  image  reconstruction,  but 
was  not  as  significant  as  improving  the  signal  strength  or  the  sampling  ratio. 

In  comparing  the  two  approaches  to  image  reconstruction,  ST Fj  is  easily  the  better  method  as 
shown  by  the  curves  remaining  more  or  less  above  unity  value  for  most  conditions  considered. 

The  advantage  of  STFi  is  greatest  for  under  =  2)  and  unity  (^  =  1)  sampling  conditions 
as  seen  by  curves  remaining  above  unity  for  much  of  the  spatial  frequency  range  out  to 

the  diffraction  limit.  Under  good  sampling  =  0.5),  STFi  approached  STFi  performance, 
particularly  at  higher  photon  counts  (N  >  200).  Becau.se  of  the  significant  bias  problem  at  the 
sampling  ratio  of  0.25,  it  was  not  determined  at  which  point  STF2  actually  outperforms  STFi 
under  very  good  sampling  conditions. 


VI 


perhaps  the  most  important  conclusion  for  this  thesis  is  that  the  SRSH  method  is  degraded  by 
a  bias.  In  tiying  to  isolate  and  characterize  an  apparent  bias  problem  with  STF  curves,  a  careful 
analysis  of  the  OTF  estimate  (.S')  showed  the  problem  of  STF  curves  exceeding  the  theoretical 
maximum  result  of  unity,  was  due  to  a  negative  bias  ((.Vr,))  in  S.  By  carefully  examining  noise 
terms  in  both  the  STF  numerator  (SS*)  and  denominator  (|.S|"),  a  useful  noise  bias  figure  of  merit 
was  developed  for  determining  non  biased  ranges  for  the  STF  curves.  The  figure  of  merit  is  the 
ratio  of  the  positive  numerator  bias  term  to  the  absolute  value  of  the  dominant  denominator  bias 
term  ( |2(s  characterize  regions  of  minimal  bias.  Basically,  the 

threshold  value  for  this  figure  of  merit  occurs  when  the  numerator  and  denominator  bias  terms  are 
identical  (value  of  0,5  for  l2(sS&Ty; ' 

5.2  Recommendations. 

The  results  from  this  thesis  clearly  point  toward  the  need  to  develop  a  non-biased  estimator. 
The  bias  problem  can  be  minimalized  by  subtracting  off  appropriately  calculated  bias  terms  at  the 
end.  However,  a  better  approach  would  be  to  develop  a  non-biased  OTF  estimate,  or  modify  the 
filtering  (Wiener  filter).  As  a  first  step,  the  negatively  biased  noise  in  S  should  be  studied  further 
to  characterize  the  influence  of  photon  count  and  sampling  ratio  values  on  STF  bias  levels. 

As  a  final  note,  some  SRSH  simulation  runs  should  be  done  using  the  minimum  variance 
control  law  for  phase  reconstruction,  since  this  was  the  control  law  used  for  the  previous  analytic 
approach  [24], 
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Appendix  A.  Selected  Numerical  Methods  From  the  Simulation. 

A.l  h'arhunen-Loeve-  Fast-Fouricr'  Transform  (h'f.-FFT)  Method  of  Phase  Screen  Generation. 

This  section  sumiriarizes  pertinent  information  from  an  unpublished  report  by  Cochran  [17]. 
The  technique  used  by  the  simulation  to  generate  random  phase  screens  involves  an  inverse  Fourier 
transform  of  filtered  complex  Gaussian  ranaom  variables,  followed  by  correction  of  low  order  spatial 
frequenies  with  a  limited  number  of  orthogonal  Karhunen  -Loeve  functions.  This  technique  is  called 
the  Karhunen-Loeve-Fast-Fourier-Transform  (KL-FFT)  Method  of  Phase  Screen  Generation.  For 
the  simulation,  the  FFT  portion  creates  a  128  x  128  pixel  array  as  represented  by  the  solid  lined 
outer  box  in  Figure  A.l.  The  x~y  axi's  are  centered  at  row  and  column  location  (65,65),  and  is 
oriented  as  shown.  KL  corrections  occur  over  a  circular  region  with  a  diameter  encompassing  129 
pixels.  A  dotted  line  on  the  lower  and  right  sides  of  the  outer  box  of  Figure  A.l  indicates  this 
circle  of  correction  extends  one  row  or  column  beyond  the  defined  boundaries  of  the  FFT  array. 
The  random  phase  screens  used  for  individual  frames  of  the  simulation  originate  from  the  center 
81  X  81  pixels  cis  indicated  by  the  inner  solid  lined  box.  Tfie  KL-FFT  method  generates  a  complex 
screen,  of  which  the  real  portion  is  used  on  odd  frames,  and  the  imaginary  portion  on  even  frames. 


A. 1.1  FFT  generated  phase  screen.  Complex  phase  screen  generation  begins  with  the 
128  X  128  pixel  array  being  loaded  with  box-muller  transform  (see  Section  A.2)  generated  complex 
gaussian  numbers,  r(  j)  =  ;i(£)-|-_/c2(*)-  The  real  and  imaginary  parts  of  c  are  zero  ntean  Gaussian 
random  variables  with  a  standard  deviation  of  one.  The  resulting  random  array  is  then  filtered  by 
the  spatial  filter  defined  by  [17:2]: 


H(f)  = 


0.1517 


if  |j|  >  0,  and 


0 


otherwise, 


(A.l) 
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Figure  A.l.  Location  of  phase  corrections. 


where  %  is  the  spatial  filter,  and  x  is  the  spatial  position.  The  FFT  estimate  of  a  coniple.\  phase 
screen  is  generated  by  the  inverse  Fourier  transform  of  this  filtered  array: 

,  (A.2) 


where  represents  the  complex  FFT  phase  screen.  The  next  step  is  to  boost  the  low  frequency 
components  of  this  screen  through  addition  of  appropriately  scaled  KL  functions. 

A. 1.2  KL  correction  to  FFT  generated  phase  screen.  Two  low  order  KL  functions  are 
used  to  correct  tilt  in  the  FFT  phase  screens.  First,  a  random  strength  is  generated  for  each  KL 
function  to  be  added  to  the  FFT  estimate  of  the  phase  screen.  The  box-muller  transform  generates 
a  gaussian  random  number,  yn  =  j/i  +jy2,  with  zero  mean  and  variance  determined  by  the  Zernike 
covariance  matrix  method,  for  scaling  the  ith  KL  function.  The  amount  of  correction  required  for 
the  ith  KL  function  is  found  by  subtracting  the  normalized  inner  product  of  the  KL  function  with 
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the  phase  screen  from  this  complex  number: 


Ex., 

Ex.,  Ej,., 


} 


(A.3) 


where  ji  is  the  complex  scaler  strength  of  the  correction  required  for  the  ith  KL  function.  The 
correction  of  the  wave  front,  for  the  KL  functions  is  then  given  by; 


^ki-fjt(t)  = (A.4) 

I 

where  is  the  KL-FFT  complex  phase  screen. 

Next,  the  real  (odd  fum  e  numbers)  or  imaginary  (even  frame  numbers)  part  of  the  center 
81  X  81  pixels  of  is  scaled  for  turbu’  .ce  strength  through  use  of  the  Fried  parameter,  r^. 

The  resulting  real  random  phase  screen  is  given  by: 

(TJe  or  Im{^i.;-//t(-?)}),  (A.5) 

' - V - ' 

center  81  x  81  pixels 

where  Tie  signifies  the  real  and  Iin  the  imaginary  part  of  Aj:  is  the  distance  between 

pixels  in  the  aperture.  This  phase  screen,  <i>,  is  used  in  the  simulation  as  the  uncorrected  phase 
front. 


A. 2  Boi-MuUer  transform. 

The  box-muller  transform  generates  two  gaussian  random  numbers  of  zero  mean  and  standard 
deviation  tr.  The  transform  is  defined  as: 


ii(-c) 

‘2(-?) 


(A.6) 

(A.7) 
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when'  (,‘i  and  (,'j  nri’  raiuloiii  nnniln'rs  Ix'lwri’n  z<to  a’lnl  one,  and  and  c-j  am  tlin  zi'ro  nn'an 
ganssian  random  nnmlnTs. 


.'1..?  h’aihiin(i)-Li>i  It  (l\l.)  Iti  prtsi  niation  of  a  Random  Phase  Scrttn. 

Tlio  Kl,  rorn-rtion  to  lli«'  ri'T  phase  screen  originates  from  tlie  techniipK'  of  using  a  seri<'s 
repres<'ntatioii  for  a  random  pliase  screen.  The  expression  for  the  random  plia.se  front  in  a  circular 
aperture  is  given  as  [17:11]: 

'i/ki{jr)  =  ^!l,l<,[r).  (A. 8) 

I 

where  'I'j./  is  Kh  generated  pha.se  screen,  A',-  is  the  ith  orthogonal  KL  function,  and  i/,  is  the 
expansion  coefficient  that  contains  the  turbulence  statistics. 

A.  i.l  ttttsis  fiinriioiis  for  a  KL  erpansion.  This  .section  summarizes  pertinent  informat  ion 
from  an  unpublished  report  by  Cochran  [17].  In  order  to  correctly  construct  a  random  phase 
sert'en.  the  basis  functions  used  must  be  ortliogonal.  One  possible  set  of  basis  functions  for  use  in 
constructing  random  pha.se  aberrations  in  a  pupil  are  the  so-called  zero-mean  normalized  Zernike 
polynomials.  These  functioils  are  defined  by: 

cosfwTr).  p  odd,  and  in  /  0, 

\/iJ(n -I- l)7?JJ’(x)sin  (mir),  p  even,  and  m  ^  0.  (A. 9) 

v/JTTTTJ”,  m  =  0. 


where  !  signifies  a  factorial  operation,  in  and  n  are  indices  of  the  Zernike  polynomials,  p  is  the 
identification  number  of  the  polynomial  according  to  .Noll's  system  [2.'):2()7  208],  and  V"  is  defined 
by: 


7?:: 


f  [(-!)»(» -i)!(kir-^’l  ] 


(A.IO) 
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where  iti  <  n,  and  n  — |ni|  is  evt'ii.  Tlie  |)rohlem  willi  Zerniko  polynomials  is  they  are  not  orthogonal 
fnnetions,  with  non-zero  correlations  existing  hetwei'ii  various  polynomials.  So  although  they  are 
ideal  for  descrihing  phase  aberrations,  they  cannot  be  used  as  basis  functions  for  a  random  ph:i.se 
screen. 

The  KL  teclmiqne  uses  an  expansion  of  the  Zernike  polynomials  to  form  orthogonal  KL 
functions.  Kach  KL  function  contains  corr<'lated  Z<'rnike  polynomials.  The  technique  involves 
first  finding  the  eigenvectors  and  eigenvalues  of  the  covariance  matrix  for  a  basis  set  of  Zernike 
polynomials.  A  KL  function  is  coti.sirncted  from  a  single  eigetivector,  whose  elements  are  u.sed 
to  scale  the  approjiriate  Zertiiki'  polynomials.  'I'he  eigenvalue  delertnines  the  variance  of  the  KL 
function. 

The  Zerniko  expansion  for  the  (th  KI.  function  is  determined  from  the  normalized  covariance 
matrix  as  follows: 

CVi=,Tfe,,  (A. 11) 


where  C  is  the  nortiialized  Zernike  covariance  matrix  for  the  basis  .set  of  zernike  polynomials,  e,-  is 
the  I'th  eigenvector  of  C,  and  <r,-  is  the  »th  eigenvalue.  The  elements  of  the  covariance  matrix  are 
calctdated  as  follows  [26:117.')],  [17:12]: 


fo.ir);i:}7(-i)'''+'-'-'-^’”V(«  +  i)(n'-bi)r(^)r(’.i±^))  , 


(A. 12) 


where  Cppi  is  the  covariance  between  numbered  Zernike  polynomials  ;i  and  p'  [25:207-208],  n  and 
m  are  indices  used  for  calculating  the  Zernike  polynomials,  r(  )  is  the  gamma  function  defined  by 


[27:380]: 


<*xp(— J’)x"“'dj-, 


(A. 13) 
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and  flaw  and  fi^a'  ar*'  Kronrfki'r  di'lla  rnnriions  d<‘(iiu-<!  .xs  [28:15;}]: 

il  irn  =  a' 

(A. 14) 

0  otherwise. 

The  /th  KL  function  is  constructed  from  the  /th  eigenvector  c,  l)y  the  following  expansion: 

=  (A.15) 

p 

where  p  is  tlie  designator  for  a  particular  Zernike  polynomial  and  c'ip  is  the  pth  comi>onent 
of  the  eigenvector  e,-.  The  phase  screen  is  found  by  summing  over  all  appropriately  scaled  KL 
functions  considered: 

=  (A.16) 

where  D  is  the  diameter  of  the  pupil,  and  pi  is  a  zero  mean  Gaussian  random  variable  with  variance 
(T?  (eigenvalue  for  »th  KL  function).  For  the  simulation,  the  two  lowest  order  KL  functions  and 
their  associated  variances  (eigenvalues)  are  calculated  as  part  of  the  setup  portion  of  the  simulation. 
The  correction  is  applied  on  odd  iterations  since  screens  are  created  two  at  a  time. 
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